Nonlinear Couplings of R-modes: Energy Transfer and Saturation Amplitudes at 

Realistic Timescales 
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Non-linear interactions among the inertial modes of a rotating fluid can be described by a network 
of coupled oscillators. We use such a description for an incompressible fluid to study the development 
of the r-mode instability of rotating neutron stars. A previous hydrodynamical simulation of the 
r-mode reported the catastrophic decay of large amplitude r-modes. We explain the dynamics and 
timescale of this decay analytically by means of a single three mode coupling. We argue that 
at realistic driving and damping rates such large amplitudes will never actually be reached. By 
numerically integrating a network of nearly 5000 coupled modes, we find that the linear growth of 
the r-mode ceases before it reaches an amplitude of around 10~^. The lowest parametric instability 
thresholds for the r-mode are calculated and it is found that the r-mode becomes unstable to modes 
with 13 < n < 15 if modes up to n = 30 are included. Using the network of coupled oscillators, 
integration times of 10® rotational periods are attainable for realistic values of driving and damping 
rates. Complicated dynamics of the modal amplitudes are observed. The initial development is 
governed by the three mode coupling with the lowest parametric instability. Subsequently a large 
number of modes are excited, which greatly decreases the linear growth rate of the r-mode. 

PACS numbers: 04.40.Dg, 04.30.Db, 97.10.Sj, 97.60. Jd 



Large scale r-modes in neutron stars become unstable 
as a result of the emission of gravitational radiation ili| 
[2|. For neutron star rotation frequencies ~ 100— lOOOHz, 
these modes emit radiation in a frequency band accessi- 
ble to gravitational wave observatories such as LIGO. 
The maximum amplitude that can be reached by these 
unstable modes not only will determine their potential 
detect ability, but also will be important for setting an 
upper limit on the rotation frequencies of neutron stars. 

The slow growth rate [3 H of the instability compared 
to the stellar rotational period hampers hydrodynamical 
simulation of the saturation of the growing r-mode am- 
plitude. Gressman et al. 3 estimate that for a 129"^ 
grid-point simulation it would take about 10^ hours of 
computational time to follow the linear growth of the 
r-mode amplitude to order unity. To circumvent this dif- 
ficulty, hydrodynamical simulations of the r-mode insta- 
bility performed to date |y| Q P have either amplified 
the radiation reaction force artificially or initialized the 
r-mode to large amplitudes. 

Schenk et al. [a] and Arras et al. Q began a pro- 
gram of using weakly nonlinear perturbation theory to 
determine the saturation amplitude of the instability. In 
this approach we represent the system as a collection of 
interacting inertial modes, and keep only second order 
interactions. Neglecting higher order interactions pre- 
supposes a small saturation amplitude, an assumption 
that can be assessed at the end. Arras et al. 9] ar- 
gued that saturation could occur at small amplitudes, 
although the argument was not as rigorous as the ap- 
proach whose first result is reported here. The advantage 
of the perturbative approach is that it permits us to focus 
on the long term evolution of mode amplitudes since we 
can analytically remove the rapid oscillations that must 
be followed in direct hydrodynamical simulations. More- 
over, our modal expansions permit considerable spatial 



resolution, which, as suggested by Arras et al. ]9||, turns 
out to be important because small scale modes affect the 
evolution substantially. 

Our starting point is Newtonian hydrodynamics for an 
inviscid self-gravitating fluid: 



dtp + V- [pv) = 
dtv + V ■ Vv = V 
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where p is the density, v is the fluid velocity, p is the 
pressure and $ is the gravitational potential. Equations 
(P) to lO are perturbed to second order to yield the am- 
plitude equations 
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where K-jg(j is the three mode coupling coefficient Q 
and eA is the rotating-frame-energy of mode A at unit 
amplitude. By specializing to a uniform density, incom- 
pressible star, we can compute Kabc analytically. We 
shall present the detailed computation of coup ling coeffi- 
cients of an incompressible star elsewhere |l3 • Although 
we have developed methods for computing k-^^q for ar- 
bitrarily fast rotation, the eigenfunctions and hence cou- 
plings and dissipation rates used here are computed to 
lowest order in rotation frequency. 

In equation Q we have augmented the inviscid hy- 
drodynamics of equations (^ to 10) by adding an imagi- 
nary part ^A to the linear modal evolution; ^a represents 
the radiation reaction force and viscous damping, which 
cause modal growth and decay respectively. The analytic 
estimates for these factors and a detailed discussion of 
their properties will be given in |10J | ; our numerical eval- 
uations of 7A agree with those obtained by Q and Q . 



By studying certain aspects of the nonlinear dynam- 
ics of the oscillator network we gain valuable insight into 
some of the challenges faced by numerical simulation. 
The problem we are considering is characterized by a 
wide range of time scales. The shortest timescale is in- 
troduced by the pulsation modes of a liquid sphere |ll| 
with frequencies 
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where the number n labels the degree of the Legendre 
functions used to describe the modal eigenfunction. It is 
this timescale that sets the time step of any direct nu- 
merical integration of the fluid equations. Notice that 
for large n, LOpui oc ^/n. Thus, for simulations with 
high spatial resolution, which we shall argue are required 
to treat the r-mode instability accurately, the time step 
must necessarily be very small, further increasing inte- 
gration times. 

The generalized r-modes or inertial modes [l^l all have 
frequencies that lie in the range 
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where 17 is the stellar angular velocity. Characteristically 
these frequencies imply timescales that are an order of 
magnitude or more slower than the pulsation timescales. 
This large mismatch in frequency will cause the inertial 
modes to interact only very weakly with the pulsation 
modes. In our investigations the pulsation modes are 
ignored. 

A third set of amplitude dependent frequencies arises 
as a result of the nonlinear couplings among the eigen- 
modes. It is this timescale that we believe both qualita- 
tively and quantitatively explains catastrophic decay of 
the r-mode observed by Gressman et al. P- 

The simplest of the class of systems described by equa- 
tion Q consists of an undamped Hamiltonian system 
that retains only one triplet of coupled modes and one 
coupling. This problem can be solved analytically in 
terms of elliptic integrals [Ijl- Although rather compli- 
cated, the solution for the modal amplitudes is periodic 
and has the generic form shown in Figure ^ The period 
depends on the initial conditions, coupling coefficient k 
and frequency detuning, Sw — Wa — wp — Wj of the in- 
teracting triplet of modes. 

For initial conditions with a large parent and two small 
daughter amplitudes, there exists a minimum parent am- 
plitude Ca = ISw/^Ky^Wjivpl below which no oscillations 
occur. Above the threshold, the period decreases roughly 
logarithmically as the parent amplitude increases. This 
threshold is the no-damping limit of the parametric in- 
stability threshold 
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discussed in detail by [Tj| [T^ O ■ ^^ ^ three mode sys- 
tem with one driven parent and two damped daughter 
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FIG. 1: (Color Online) Generic solution of the three mode sys- 
tem starting with a large initial amplitude in the r-mode and 
negligible amplitudes in the two coupling daughter modes. 



modes, this quantity sets the minimum amplitude the 
parent mode has to reach before the daughters are ex- 
cited. 

For systems with substantial initial r-mode amplitudes, 
the time that elapses before the r-mode first declines 
rapidly via nonlinear three mode effects depends on the 
amplitude itself. The evolution of the r-modes in a sys- 
tem of three coupled modes starting at various initial am- 
plitudes is shown in Figure|21 The coupling used to make 
these plots was selected to yield the minimum paramet- 
ric instability threshold possible with the resolution used 
in the simulations by Gressman et al. Q- We estimated 
this resolution by considering that cartesian grids with 
129'^ and 65'^ grid points will be able to resolve modes 
with ri = 15 by having at least 4 points per wavelength 
even if these modes will not be represented accurately. 

Note that the normalization convention chosen 
in jal results in an energy per unit amplitude of 
E = ^a^R^Q^p which differs from the normalization of 
our amplitudes Cq to £■ = ^AIR^fl^ at unit amplitude. 
The conversion between the two amplitudes is given by 
Ca — ■\/3/567ra w O.lSa . In order to get the time axis in 
our figures in ms to compare with the results of [5| you 
need to multiply by a factor of 1/217 = 0.099ms. 

The three mode system was started with initial ampli- 
tudes for the r-mode that had energies equivalent to those 
used to generate Figure 3 left of Q • The daughter modes 
were assumed to have initial amplitudes of r-^ 10~^. 

The top four results of Figure |21 correspond to the am- 
plitudes plotted in the left part of Figure 3 of la] and 
are qualitatively similar to what was found in pj. In 
5] these amplitudes result in a difference between de- 
cay times of At ^ 30ms for the "pumped up" state; our 
three mode coupling example results in decay times dif- 
fering by Ai ~ 28ms for the same initial amplitudes. Dif- 
ferences in the initial amplitudes of the daughter modes 
and detuning conditions, resulting from different eigen- 
f unctions in the two problems, may change the numer- 
ical result slightly. However, considering that our cal- 
culation is based on a rough knowledge of the resolu- 
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FIG. 2: Hamiltonian three mode system with large initial am- 
plitude. Comparing the times to the first decay of the parent 
mode as a function of initial amplitude. Graphs computed 
for the coupling with the lowest parametric threshold with 
n < 30, the coupling details are k = 0.191, Sw = —1.91 x 10~^ 
the frequencies of the modes in the co-rotating frame are 
{-0.219985, 1/3, -0.113349}2n 



tion of the simulation grid used in Q , and employs the 
coupling coefficients and detuning of the inertial modes 
of a nearly spherical incompressible star with negligible 
damping and driving, the simple three mode model cor- 
responds to the hydrodynamical simulations remarkably 
well. This indicates that what was observed in [3 is 
the decay of an unphysically "pumped up" state to just 
a few accessible modes with initially small amplitudes. 
It should be very interesting to compare the ratios of 
the low frequency modes visible in the Fourier spectra 
of Gressman et al.'s decaying solution on the right hand 
plot of their Figure 3 to those used to produce Figure |21 
Unfortunately the simulation |3| was stopped before inte- 
grating long enough to determine these lowest frequency 
components accurately. 

Therefore the simulations reported in lq| evidently cor- 
respond to a nonlinear system involving only a few im- 
portant modal couplings to the r-mode. In actuality the 
evolution of the r-mode instability from an initially small 
amplitude may be much more complicated. 

By calculating the parametric instability thresholds of 
the r-mode for its 146998 direct couphngs to modes with 
n < 30 we find that the lowest two parametric insta- 
bility thresholds are Ca = 1.6 x 10^^ for the coupling 
to modes (n — 13, m — 3) and {n — 14, m — 1) and 
Ca = 1.2 X 10~^ for the coupling to modes with (n = 15, 
m = 7) and (n = 14, m — 5), where m is the azimuthal 
quantum number of the mode. In our simulations we 
find that both thresholds are significant in the evolution 
of the r-mode. Each time the mode amplitude passes a 
certain threshold there is a dramatic change in its dy- 
namics, generally decreasing its growth rate. A detailed 
discussion of these simulations and the exploration of the 
effect of different damping and driving rates will be pre- 
sented in a separate paper. An instructive example is 
shown in Figure |2| 
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FIG. 3: (Color Online) Evolution of the mode amplitudes 
for a coupling network with 4959 modes with n upto 30, the 
1306999 couplings with Sw < 0.002 were included, n = m + 1 
r-modes are shown in black. The mode that reaches the 
largest amplitude is the most unstable r-mode. The lowest 
two parametric instability thresholds for this mode are in- 
dicated as horizontal dashed lines. The total energy of the 
system is shown above. 



The gravitational driving for this simulation corre- 
sponds to the radiation reaction force experienced by 
a star rotating at f2 = O.Sly/nGp . Damping is due 
to shear viscosity calculated for a star with temperature 
T = 5 X 10'^ K which implies that any mode with n > 25 is 
capable of individually damping out the r-mode growth. 

Following an initial period of linear growth, the r-mode 
reaches its lowest parametric instability threshold in our 
network. The dynamics are then governed by the this 
particular unstable three mode coupling. In Figure El 
the first three bumps in the r-mode amplitude are a key 
signature of a three mode coupling in which the two 
daughter modes do not damp strongly enough to stop 
the growth of the parent mode amplitude. After the 
second parametric instability threshold is reached, many 
new pathways to distribute energy become available, and 
there is no longer a predictive way of describing the evo- 
lution of the system. We do observe that a number of 
modes with n > 24 are excited at this stage. These modes 
result in considerable dissipation that rapidly damp the 
r-mode. 

During the simulation we reached an integration time 
of ^ 10^/217. Existing hydrodynamical simulations last 
about 10^/2r2, which is not long enough to follow the en- 
ergy transfer from the r-mode to other fluid modes start- 
ing from small amplitudes . The time scales involved 



for such a transfer require integration times of at least 
10^/217. For a similar reason greatly amplifying the ra- 
diation reaction force will also inflate the r-mode unnat- 
urally, before it has time to react with other fluid modes, 
so it may attain unphysically large amplitudes. Because 
of our long integration timescales, we are able to follow 
modal interactions and decay with realistic radiation re- 
action forces. 

Although the linear growth of the r-mode ceases after it 
has passed the first two parametric instabilities, there is 
no rigorous argument that the r-mode cannot continue to 
grow, and slowly pump energy into the network of modes 
until establishing a cascade as suggested by Arras et al. 
[3|. This growth rate will necessarily be slower than the 
linear growth rate and require longer integration times. 
If longer integration times can be attained, it may be 
possible to check if a steady state cascade proposed by 
Arras et al. j^ is realized in the limit of small damping. 
However even our simulation may involve too few modes 
to check this effectively. 

The resolution used in simulating an oscillator net may 
influence its dynamics. This is best illustrated by means 
of an example. During our first investigations we used 
a network with n < 12. We found that our integrations 
diverged soon after the r-mode decayed for the first time 
regardless of how strongly the daughter modes damped. 
Energy "bounced" back off the bottom end of the spec- 
trum before the solution diverged, indicating that too few 
modes had been included and/or that our small ampli- 
tude assumption failed. 

The reason for this failure is that the lowest parametric 
instability threshold for the subset of modes with n < 12 
for our model is Cq > 0.0058. Once this r-mode ampli- 
tude was reached via means of linear growth, the three 
mode coupling time scales were much faster than the 
modal damping rates, causing energy to be transferred 



rapidly to larger n modes, and resulting in the diver- 
gence. 

In a simulation with coarse resolution, the possible 
ways in which the r-mode can transfer energy to other 
modes is limited. The altered dynamics of such a trun- 
cated system tends to promote unphysically large ampli- 
tudes. Our much larger network of n < 30 did not exhibit 
this pathology. 

To summarize, we have shown that simulations start- 
ing from small amplitudes in all modes (except possibly 
the r-mode) will initially be governed by the dynamics 
of three mode coupling. We demonstrate that the rapid 
decay of the r-mode observed in numerical simulations 
@| is consistent with that predicted from perturbation 
theory in a truncated three mode system. 

However, we further argue that if realistic values of 
the radiation reaction driving are considered and a large 
network of interacting modes is followed, the r-mode be- 
gins to interact nonlinearly with other fluid modes at an 
amplitude of about 10~^, resulting in much longer decay 
timescales and smaller amplitudes than those attained in 
existing hydrodynamical simulations. We also find that 
if the damping is weak enough, the r-mode couples to a 
large number of modes, suggesting that a possible cas- 
cade may be established, although our simulations cov- 
ering r^ 10^ stellar rotation periods were evidently too 
short to reach such a state. This will be the subject of 
further investigation. 
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